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We consider the astrophysical consequences of the self interacting dark matter 
(SIDM) scenario for a general velocity dependent cross section per unit mass 
which varies as some power of velocity: a dm — °~o ( v /vo)~ a . Accretion of SIDM 
+2 • onto seed black holes can produce supermassive black holes that are too large 

for certain combinations of o~qVq an d a, a fact which is used to obtain a new 
constraint on the dark matter interaction. Constraints due to other astrophysical 
considerations are presented and previous constraints for a constant cross section 
are generalized. The black hole constraint is extremely sensitive to the cusp slope 
a, of the inner density profile p ~ r~ a of dark halos. For the most probable value 
of a = 1.3, we find that there exists a tiny region in the parameter space for SIDM 
properties, with a rj 0.5 and (cr /l cm 2 g _1 ) (t> /100 km s _1 ) a ~ 0.5, such that 
all constraints are satisfied. However, the adiabatic compression of the dark halo 
by baryons as they cool and contract in normal galaxies yields a steeper cusp, p ~ 
r~ a ' . We find that in both the highly collisional and collisionless limits, invariance 
arguments require a' = a+ f_~ , where a and a' are the inner profile slope of 
the dark halo before and after compression, respectively. This gives the tighter 
constraint (cr /i cm 2 g -1 ) (t> /100 km s _1 ) a < 0.02, which would exclude SIDM 
as a possible solution to the purported problems with CDM on subgalactic scales 
in the absence of other dynamical processes. Nevertheless, SIDM with parameters 
consistent with this stronger constraint, can explain the ubiquity of supermassive 
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black holes in the centers of galaxies. A "best fit" model is presented with a = 
and (o" /l cm 2 g" 1 ) = 0.02, which reproduces the supermassive black hole masses 
and their observed correlations with the velocity dispersion of the host bulges. 
Specifically, the approximately fourth power dependence of black hole mass on 
galactic velocity dispersion is a direct consequence of the power spectrum of 
primeval perturbations having an index of n pa —2 and the value of a. Although 
the dark matter collision rates for this model are too small to directly remedy 
problems with CDM, mergers between dark halos harboring supermassive black 
holes at high redshift could ameliorate the cuspy halo problem. This scenario also 
explains the lack of comparable supermassive black holes in bulgeless galaxies like 
M33. 

Subject headings: Dark matter - galaxies: formation, halos - black hole physics 
-cosmology: theory 

1. Introduction 

Self interacting dark matter (SIDM) has been proposed by Spergel & Steinhardt (2000) 
to remedy purported problems with the cold dark matter (CDM) family of cosomological 
models on subgalactic scales . The apparent conflicts between numerical simulations of 
CDM and observations from dwarf galaxies to clusters seem to indicate that CDM halos 
are too centrally concentrated (see Wandelt et al. 2000 and references therein). Of these 
discrepancies, two problems are particularly significant. First, the inner density profiles of 
CDM halos diverge as p oc r~ a with a w 1.3 ± 0.2 (Navarro, Frenk, & White 1996, 1997; 
Fukushige & Makino 1997; Moore et al. 1999; Subramanian, Cen, & Ostriker 2000; Jing 
& Suto 2000; Ghigna et al. 2000; Klypin et al. 2001; Fukushige & Makino 2001a, 2001b), 
possibly conflicting with profiles deduced from rotation curve observations of dark matter 
dominated dwarf and low surface brightness galaxies (Flores & Primack 1994; Dalcanton 
& Bernstein 1997; De Blok & McGaugh 1997; Swaters, Madore & Trewhella 2000; though 
see van den Bosch et al. 2000 and van den Bosch & Swaters 2001). Second, numerical 
simulations of the CDM scenario indicate an excess number of small scale structures when 
compared with the number of dwarf galaxies in the Local Group or satellites in galactic 
halos (Moore et al. 1999; Klypin et al. 1999). Endowing dark matter with a finite cross 
section for elastic scattering allows heat transfer to occur over a Hubble time, resulting in less 
concentrated structures, flattened density profiles, and fewer substructure satellites (Burkert 
2000; Yoshida et al. 2000b; Dave et al. 2000; Wandelt et al. 2000). 

If a dm = is the cross section per unit mass of a dark matter particle, then the 
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mean free path is A = 1/p<Jdm, where p is the density of dark matter. An optical depth 
can be defined r = r/A, which distinguishes different physical regimes. In regions which are 
optically thick t > 1, the dark matter behaves as a fluid; whereas, for regions which are 
optically thin r < 1, the dynamics closely resembles two body relaxation in globular clusters. 
The range of cross sections consistent with experimental constraints, yet giving the required 
rates of evolution, imply that dark halos will probably be optically thin r (r s ) < 1 at their 
characteristic scale radius r s (Wandelt et al. 2000; Dave et al. 2000), but of course can be 
optically thick in the inner regions r<r s where the densities are higher (provided a > 1). 

Numerical investigations of SIDM halo evolution have been carried out by several groups. 
Yoshida et al. (2000a) and Moore et al. (2000) simulated SIDM halos in the optically thick 
or fluid limit, which, as mentioned above, is probably not the relevant scenario for SIDM. 
Burkert (2000) and Kochanek & White (2000) simulated isolated halos employing a range of 
cross sections nearly consistent with the optically thin requirement. Both groups found that 
halos develop shallow cores for a modest period of time before the onset of core collapse, 
but disagree on the core collapse timescale. Dave et al. (2000) and Yoshida et al. (2000b) 
simulated optically thin halos in cosmological settings, which include infall and merging. 
They observed evolution towards reduced central densities and shallower inner profiles, effects 
increasing with increasing cross section; and neither saw any evidence for core collapse. Their 
results are in broad agreement with each other and span the range from dwarf galaxies to 
clusters. Alternative to N-body techniques, Hannestead (2000) and Firmani, D'Onghia, 
& Chincarini (2001) have carried out integrations of the collisional Boltzmann equation. 
Firmani et al.'s simulation differs from previous investigations in that they considered a 
velocity dependent cross section a oc i> _1 . Also Balberg, Shapiro & Inagaki (2001) examined 
SIDM halo evolution via a time-dependent gravothermal numerical calculation. All three 
find the development of shallower central slopes and less concentrated cores than CDM would 
produce, consistent with the quoted N-body results, although Balberg et al. (2001) obtained 
a collapse timescale roughly an order of magnitude larger than that seen by Burkert (2000) 
and Kochanek & White (2000). In sum, the most relevant numerical studies performed 
to date indicate that, for appropriate dark matter scattering cross sections, collisions can 
effectively reduce dark matter central densities for normal galaxies. 

A number of authors have obtained constraints on <tdm from analytical and semi- 
analytical arguments. Strong lensing events are extremely sensitive to the inner profiles and 
shapes of dark halos in clusters of galaxies, and are thus a powerful probe of SIDM, which 
has enabled several authors to place constraints on a DM (Miralda-Escude 2000; Wyithe, 
Turner & Spergel 2001; Meneghetti et al. 2001). Mo and Mao (2000) determined what value 
of (Jdm would produce a correct Tully-Fisher relation for SIDM halos. Further constraints 
can be placed on SIDM because of the existence of subhalos in larger halos: heat transfer to 
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a cool subhalo from the hotter halo could evaporate subhalos conflicting with observations if 
the evaporation were too efficient. Gnedin and Ostriker (2001) constrained a DM by requir- 
ing that the dark halos of elliptical galaxies in clusters survive until the present. Finally, 
Ostriker (2000) showed that accretion of SIDM onto seed black holes could produce super- 
massive black holes in the observed range and with the observed mass scaling Mbh oc v 4-5 
(Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000; Merritt & Ferrarese 
2001), and he placed an upper limit on cr DM so as to avoid the formation of central black 
holes that are too large. 

In this paper we consider the astrophysical consequences of a general velocity dependent 
cross section, as has been discussed by various authors (Yoshida et al. 2000b; Dave et al. 
2000; Hogan & Dalcanton 2000; Gnedin & Ostriker 2001; Firmani et al. 2001a,2001b; Wyithe 
et al. 2001; Hui 2001). We refine the black hole constraint of Ostriker (2000), present new 
constraints due to other astrophysical considerations, and generalize some of the constraints 
discussed above for a non-constant cross section. The ability of SIDM to reproduce the 
observed properties of the galactic supermassive black hole population is also explored. In 
§2 we determine the black hole mass that will be grown in a given dark halo. We discuss in 
§3 the sensitivity of this result to inner profile flattening. We consider in §4 the scaling of 
black hole mass with halo velocity dispersion. In §5 we obtain upper limits on the strength 
of the dark matter interaction from the observed masses of supermassive black holes. The 
effect of the merger history of dark halos is discussed in §6. New constraints on the dark 
matter interaction are presented in §7 and we then summarize and generalize previous ones. 
We conclude in §8. Throughout this work we adopt a ACDM cosmology with f2 = 0.3, 
tt b = 0.045, tt A = 0.7, h 100 = 0.65, a 8 = 0.9 (see e.g. Bahcall et al. 1999; Eisentstein & Hu 
1999) 



2. SIDM Black Holes 

A fundamental assumption of our black hole growth scenario is that SIDM dark halos 
have singular density profiles p ~ r~ a with 1 < a < 2 immediately after collapse. Before 
discussing black hole formation, we digress briefly to justify this assumption. The essential 
point is that the collapse occurs on a dynamical time scale (~ 10 7 yr for galactic halos), 
whereas the subsequent flattening of the inner profile due to dark matter collisions occurs as 
the halo evolves secularly on a much longer relaxation timescale (~ 10 10 yr). The singular 
density profiles found in numerical simulations of CDM can been understood with the help of 
two different analytical pictures. The first is the spherical secondary infall paradigm (Gunn & 
Gott 1972; Fillmore & Goldreich 1985; Bertschinger 1985; Hoffman & Shaham 1985; Hoffman 
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1988; Subramanian et al. 1999), whereby a perturbation from the homogeneous cosmological 
background collapses from the inside out, with subsequent shells of matter collapsing at 
succesively later times. Bertschinger (1985) demonstrated that the same singular density 
profile is to be expected for both collisionless and collisional gases, for the inner parts of the 
self-similar solutions of the spherical collapse problem. The alternative picture for explaining 
cuspy profiles is violent relaxation (Lynden-Bell 1967), the conjecture being that particles 
will transfer energy to one another as they move through the rapidly varying potential of the 
collapsing system, resulting in a singular isothermal density profile p oc r~ 2 . Whatever the 
dominant mechanism, both spherical collapse and violent relaxation occur on the dynamical 
timescale of the halo, which will be much shorter than the timescale for heat transport (via 
dark matter collisions) which flattens the inner halo (see §3). Kochanek and White (2000) 
tested this hypothesis by simulating the formation of SIDM dark halos. They compared 
the density profiles resulting from the collapse of a top-hat overdensity with and without 
collisions and found that the final equilibrium profiles are similar and can be fit by cuspy 
profiles in both cases. Both the analytical frameworks and Kochanek and White's numerical 
experiment make it highly unlikely that SIDM dark halo profiles are flat at early times. 

To this end, we assume the density profile at early times will be a generalized Navarro- 
Frenk- White (NFW) (Navarro et al. 1996,1997) or Zhao profile (Zhao 1996) 

p( r )=Ps(Pj r. ~ ) Y - a , (2-1) 
\r s J (1 + r/r,) 

so that in the inner regions r<r s the density is given by p (r) w p s (jrj with 1 < a < 2 
but most likely a — 1.3 ± 0.2 (Navarro et al. 1996, 1997; Fukushige & Makino 1997; Moore 
et al. 1999; Subramanian et al. 2000; Jing & Suto 2000; Ghigna et al. 2000; Klypin et al. 
2001; Fukushige & Makino 2001a, 2001b). Imposing hydrostatic equilibrium for this profile 
gives the velocity dispersion 

/ r \ 2-a 

v 2 {r)=v 2 a l-J ; r<r s , (2-2a) 

where 
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v 2 s = pGp s rj ; » = - (2-2b) 

(6 — a) (a — Ij 



Recent high resolution numerical studies of pre-galactic objects indicate that massive 
stars (~ 100 M Q ) form at the centers of progenitor dark matter halos at redshifts z— 15 — 20 
(Abel, Bryan, & Norman 2000, 2001). Subsequent evolution will end in supernovae on Myr 
timescales, leaving a population of remnant (~ 50 M ) seed black holes, long before the 
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first galactic mass dark halos have formed (Madau & Rees 2001). Therefore, the material 
that accretes and merges to form galactic mass dark halos will be well seeded with black 
holes already accreting dark matter on dynamical timescales shorter than the timescale 
for SIDM collisional evolution. Following Ostriker (2000), we consider the quasi-spherical 
accretion of dark matter onto a single seed black hole at the center of a galactic dark matter 
halo. Clearly, the assumption of a single seed black hole at the halo center is an over- 
simplification. However, given the rapid phase of initial growth (see below), one black hole 
is likely to dominate and will eat or eject the others, so that this complication should not 
alter our estimates of the final black hole mass. 

The dark matter is treated as an adiabatic gas, which is valid provided the optical depth 
diverges for r << r s (an assumption we check below). The density and velocity dispersion 
can be written 

,(r)»H*)T r * r ' , P-3) 

and 

C Y( 2 _ a ~ , (2-4) 

v 2 a f r c < r < r s 



(2-5) 



where r c is the "capture" or Bondi accretion redius 

GM BH 



v 2 (r c ) ' 



The quantities p c and v c can be expressed in terms of p s and v s by requiring p and v be 
continuous. 

Consider a cross section per unit mass u DM , which varies as some power of velocity 

O-DM (v) = (To ( — J , (2-6) 

where <7o, Vq, and a are determined by the fundamental physics of the interaction. Then the 
optical depth is 

f f r \ {a - 1)/2 
I t.. I — I r <~ r .. 

(2-7a) 




where 

= PcTcOq (v c /v y a , (2- 7b) 
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and 

P = a(l-a/2) + a. (2-7c) 

For our halo to be optically thick in the inner regions (r 3> 1 for r <C r c and r <C r s ), we 
must have a < 1 and (3 > 1. The latter condition is satisfied for all a > 0, provided that 
1 < a < 2 ; however, the former implies a hydrodynamic treatment of the accretion flow is 
only valid provided the cross section velocity dependence is not too steep. Hence, the black 
hole mass we derive below only applies to cross sections with < a < 1. 

The characteristic size of our accreting hydrodynamic system during this period of rapid 
growth is r c (eqn. 2-5). The black hole will grow via Bondi-like accretion until the accretion 
radius is comparable to the mean free path of the dark matter, r (r c ) = 1 . Then a transition 
will occur to a much slower, diffusively limited growth. This allows us to obtain a rough 
estimate of the black hole mass grown during the optically thick phase 

M BH « ^ Ps r 3 s rh . (2-8) 

Thus, the black hole mass will be of order the halo mass, p s r 3 , times the optical depth 
at the characteristic scale, t s , raised to a power which depends on the inner density profile 
and the velocity dependence of the scattering. Figure 1 shows the dependence of the black 
hole mass on the inner profile slope for a Milky Way size dark halo, M^aio = 3.0 x 10 12 M , 
and several different values of the exponent, a, of the cross section velocity dependence. The 
black hole mass depends on the combination uqVq of the cross section parameters in eqn. 
(2-6). In Figure 1 we set o\v f Q0 = 1 where o\ = (<7o/l cm - 2 S" 1 ) an d ^ioo = (i>o/100 km s _1 ). 
The characteristic density, radius, and velocity dispersion, p s , r s , and v s have been calculated 
from the dark halo concentrations using the routine made publicly available by Eke, Navarro, 
& Steinmetz (2001) (hereafter ENS), and described in Appendix A. 

Once the accretion radius extends into the optically thin region of the dark halo the fluid 
approximation is no longer valid. The subsequent slow growth of the black hole proceeds 
as dark matter particles are scattered into the loss cone. This phase of growth was treated 
by Ostriker (2000) for a velocity independent cross section. Loss cone accretion can grow 
black holes comparable to eqn. (2-8) (for a general velocity dependent a dm)- However, this 
assumes a cuspy profile extends into the innermost regions of the halo for a Hubble time, 
which will not be the case if halos are significantly flattened by heat transfer during that 
same time interval. Conversely, accretion from the optically thick region of the halo grows 
black holes nearly instantaneously in comparison to cosmological timescales, as we will see 
in the next section. Since we have neglected the optically thin phase of growth and also any 
contribution to the mass from baryons, eqn. (2-8) should be regarded as a lower limit when 
compared to observations. Finally, it should be noted that the black hole mass calculated 
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here is that grown from a single dark halo. We will consider the effect of the dark halo 
merger history on this estimate in §6. 



3. Inner Profile Flattening 

From Figure 1 it is apparent that the black hole mass in eqn. (2-8) may be very sensitive 
to the value of the inner profile exponent a. It is thus prudent to determine how physical 
processes not considered above might alter the inner profile. In what follows we consider the 
effect of the accretion flow and heat transfer. 

If, to lowest order, we take the accretion to be spherical, then, in the inner regions 
r <^ r c <^ r s , gravity dominates over pressure and the particles are in free fall with the radial 
velocity u oc r~ 1//2 . Then M = const, implies p oc r~ 3 / 2 . This density profile will interpolate 
smoothly with the most probable value of the inner profile slope a — 1.3 ± 0.2 at r ~ r c . 

Now consider the effect of heat transfer. We have just seen that u oc r -1 / 2 or T oc r _1 
for r <C r c , so that the inner regions of the flow are dynamically hot. Also the density profile 
in eqn. (2-1) produces a temperature inversion, so that heat will flow from the outer halo 
inward flattening the inner profile. We are thus led to consider how heat flow outward from 
the accretion flow r<r c and inward from the outer halo r S> r c alters the density profile 
near the temperature minimum at r ~ r c . Taking note of the fact that the transport behavior 
differs under optically thick and thin conditions, we consider these cases separately. Figure 
2, presents a cartoon to illustrate the density and temperature profiles and the direction of 
heat flow for an SIDM dark halo with a black hole accreting from the inner regions. 

For r << r c (r ^> 1), the dark matter behaves like a fluid and hence heat transfer can 
be described by a diffusion equation 

T dS = ^d_ r 2 m K dT 
dt pr 2 dr p <9r ' 

where S is the entropy, m p is the dark matter particle mass, and T is the temperature defined 
as \ksT = \m p v 2 . The coefficient of thermal conductivity for the monatomic gas is 

p c v ( k B T \ 1/2 

k ~ c v \ — v = — , (3-2) 

m p m p \cr DM m p J 

where c v is the specific heat per particle and A is the mean free path. The entropy S relates 
to the equation of state P = -£-k B T oc p 7 . It can be written 

^dS d m d 1 , , d c v T . . 

T— = —C V T + P-r— = O 1 —-^—r. 3-3 

dt dt dtn H dtp-r- 1 v ; 
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Combining eqns. (3-1), (3-2), and (3-3) and de-dimensionalizing them using the character- 
istic size r c , density p (r c ), and temperature T c = mpV k ^ Tc ^ , we obtain for the unit of time 

^thermal = T ( r c) tdyn ( r c) • (3"4) 

Here, thermal is the heat conduction timescale, and t dyn (r) = r/v (r) is the dynamical time. 
Since the dark matter will typically be optically thick at r c , t t h erma i 3> td yn (r c ). This 
condition breaks down in the final stages of black hole growth when the accretion radius r c 
approaches the region of the halo of optical depth unity r ~ 1. Again, taking the accretion 
flow to be spherical, one can show that the accretion radius grows (as the mass of the black 
hole grows) at the local sound speed 

^ = 2(a-l)v(r c ). (3-5) 

Then,the speed at which thermal disturbances propagate is 

r c v (r c ) dr c 

1 'thermal ~ T ~ 7 x < "jT, (<3-OJ 

^thermal T vc) 

which is much slower than the growth of the accretion radius. We conclude that heat 
conduction from the the accretion flow outward cannot alter the density profile fast enough 
to halt the growth of the black hole, until the accretion radius grows into the region of the 
halo where r (r c ) ~ 1. 

Now we consider heat conduction inward from the optically thin halo. Integrating eqn. 
(3-5) gives r c (t); and we can solve for the black hole growth time tsH, which is the time for 
the accretion radius, r c , to grow until r (r c ) = 1 : 

t B H = r(r s )^t dyn (r s ). (3-7) 



y 100 ~ 



For the range of parameters which satisfy all the astrophysical constraints (see §7), <Jivf 
and < a < 1, r s < 0.02 on galactic scales. Thus, from eqn. (3-7) (tBii/tdyn (r s )) ^ 1, and 
the dark matter induced phase of black hole growth will be quite rapid, essentially on the 
timescale of the collapse of the dark matter halo. From the dwarf galaxy to cluster scale 
with a = 1.5, we find 10 5 yr < tsH ^ 10 7 yr for < a < 4. In particular, for the case a = 1 
discussed in the literature (Yoshida et al. 2000b; Dave et al. 2000; Hogan & Dalcanton 2000; 
Gnedin & Ostriker 2001; Firmani et al. 2001a,2001b; Wyithe et al. 2001; Hui 2001), t BH is 
independent of the scale of the halo and is given by 

a v n5 / a \ ( v 



tBH = ^ = 5.7 x 1Q6 yr p. . (3-8) 

jj,G J \cm 2 g' 1 J V100 km s~ 1 J K ' 
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As mentioned above, the post-collapse profile of the dark halo has a temperature inver- 
sion. Ultimately, the conduction of heat inward from the outer halo will flatten the inner 
profile and a constant density core will form. This evolution will take place on the relaxation 
timescale (Kochanek & White 2000; Burkert 2000; Quinlan 1996, Balberg et al. 2001) 

Uei (r s ) = -t dyn (r s ) = t BH r^ 1+m ^^ » t BH . (3-9) 

The relaxation time is of order a Hubble time for the optically thin halos considered, so the 
black holes are effectively grown instantaneously in comparison to the relaxation time and 
other cosmological timescales. Although SIDM serves to flatten central density profiles after 
relaxation timescales, before this occurs 10 6 — 10 8 M & black holes will have grown in galaxies 
by optically thick accretion provided their dark halos had a cuspy profile for ~ 10 6 yr of their 
history. 



4. Scaling Relations 

Figure 3 shows plots of the black hole mass from eqn. (2-8) versus dark halo cir- 
cular velocity for different values of a, where we have taken <Tiw" 00 = 1. These were 
generated using the ENS scaling relations described in Appendix A (ENS 2001). The 
dashed line is the M B h — v relation determined by Merritt & Ferrarese (2001), M BB — 
1.3 x 10 8 M (cr c /200 km s^ 1 ) 4 ' 72 , where a c is the central velocity dispersion of the bulge. 
The dotted line is the shallower M BH — v relation derived by Gebhardt et al. (2000), 
M B h = 1-2 x 10 8 M Q (<r e /200 km s^ 1 ) 3 ' 75 , where a e is the luminosity-weighted line-of-sight 
velocity dispersion within the half-light radius. Note that, for the sake of comparison to our 
model, we have extrapolated these observed M BB — v to all dark halo scales, although black 
holes have only been observed in galactic scale systems. In order to compare SIDM black 
holes to observations we must relate the circular velocity of the dark halo to these measured 
bulge velocity dispersions, a c and a e (or equivalently we could relate the mass of the halo to 
the mass of the bulge). A determination of these relationships is beyond the scope of this 
work, but for the sake of comparison we naively assume simple proportionalities, a c oc v circ 
and a e oc v circ , between the quantities. The constants are determined from the Milky Way 
for which o c = 100 km s _1 (Merritt & Ferrarese 2001) , a e = 75 km s _1 (Gebhardt et al. 
2000), and v circ = 220 km s" 1 (Binney & Tremaine 1987). 

It is apparent from Figure 3 that increasing the exponent a tends to flatten the M BH —v 
relation for black holes grown from SIDM. Ostriker (2000) considered the simplest case, 
a = 0, and found M BH oc v 4 , in satisfactory agreement with observations. Figure 3 also shows 
our "best fit" model to the magnitude and slope of the two observed M BH — v relations, 



- 11 - 



which has a = 1.74, o-iv% 00 = 0.02, and a = 0. Setting a = reproduces the slope of 
the observed relations, while the other parameters, o"if" 00 and a, have been chosen to be 
consistent with the discussion in the next section. However, it should be kept in mind that 
other combinations of a and (Jivf Q0 would yield a similar Mbh — v relation because of a 
degeneracy in these two quantities. 

The scaling in Figure 3 can be understood in a cosmological context as follows. Suppose 
a halo of mass M s collapses to form a virialized object with characteristic size r s , density 
p s , and velocity dispersion v s . For a power law power spectrum of density fluctuations 
P (k) oc k n with a spectral index —3 < n < 1, simple scaling arguments from linear theory 
predict (Peebles 1980; Padmanabhan & Subramanian 1992; Padmanabhan 1993) 

p s <xvf n+ VI {n ~ l \ rsccv; 2 ^^. (4-1) 

Use of these relations in eqn. (2-8) gives the scaling of black hole mass with velocity disper- 
sion. If a = 

4(ra+2) (3-a) 12 

Mbh oc u a lB_1) (< * _1) ^ ■ (4-2) 

On galactic scales the spectral index n ~ —2, which gives the approximate scaling Mbh oc v^, 
independent of a. For general a and taking ct = 1.5 and n = — 2 

M^oc^^ 17 , (4-3) 

which roughly explains why the M B h — v relation is flatter for larger a. 



5. Black Hole Constraints 



Thus far, we have completely neglected baryons in our discussion of supermassive black 
holes grown from accretion of dark matter. Baryons will alter the estimates of M B h in two 
ways. First, accretion of baryons will add to the mass of the central black hole so that 
the black hole mass estimated in §2 should be considered a lower limit when compared to 
observations. Second and perhaps more significant, as baryons cool and contract they will 
condense in the inner regions compressing the dark halo and resulting in a steeper inner 
density profile. If this compression occurs before the flattening caused by dark matter colli- 
sions, then the black holes grown will be significantly larger because of the steep dependence 
of M B h on a in Figure 1. Observations indicate that in general, bulgeless galaxies do not 
appear to contain massive black holes comparable with those found in ellipticals of compara- 
ble total mass (Richstone et al. 1998; Kormendy 2000; Gebhardt et al. 2001). In the black 
hole growth scenario described here, the explanation for this could be that galaxies with 
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inner regions dominated by a stellar bulge, will have compressed their central dark halos 
thus increasing a, whereas, compression will be less significant in bulgeless systems. 

With regards to constraining the dark matter interaction, the question arises as to 
whether eqn. (2-8) should be compared to bulge or bulgeless systems. In this section, we 
first obtain a clean constraint by comparing to the bulgeless system M33, where the effects 
of baryons can be neglected. In Appendix B, the effect of baryonic infall and compression 
on the inner profile slope of the dark halo is determined, which is used to obtain another 
black hole constraint by comparing to the bulge systems that lie on the Mbh — v relation. 

5.1. A Supermassive Black Hole in M33? 

M33 is a normal low luminosity dark matter dominated spiral (Scd) which lacks a 
significant bulge. Its rotation curve is dominated by its dark halo from 3 kpc outward, so 
that it is relatively easy to disentangle dark and luminous matter and obtain an estimate 
for the mass of the halo (Persic, Salucci, & Stel 1996; Corbelli & Salucci 2000). Corbelli & 
Salucci (2000) analyzed the rotation curve of M33's H I disk out to 13 disc scale lengths (16 
kpc). From their outermost data points, v C i rc as 125 km s _1 , which we take as the circular 
velocity of the dark halo, allowing us to deduce its mass, M ha i « 5.1 x 10 11 M , from the 
scaling relations in Appendix A. M33 also has a notoriously small upper limit on the mass of 
its central black hole. Observations from the Space Telescope Imaging Spectrograph (STIS) 
on the Hubble Space Telescope place an upper limit of Mbh ^ 1500 M Q (Kormendy et al. 
2001; though see Merritt et al. 2001). 

The masses of the dark halo M ha i Q and black hole M BH allow us to obtain a constraint 
on the dark matter interaction from eqn. (2-8). Specifically, an upper limit on the quantity 
<tiViq can be obtained as a function of a: 

^i<oo £ ^ „ = 23 - 9 (0.037) a , 2.34 (0.051) a , 0.13 (0.10) a 

p s r s (1 cm* g-i) 

(5-1) 

where p s , r s , and v s are determined by M ha [ from the scaling relations, and we have set 
a = 1.1, 1.3, and 1.5 respectively to obtain the last equality. The accretion of baryons 
will increase the black hole mass above our estimate in eqn. (2-8) so that eqn. (5-1) is an 
overestimate and the constraint could be tighter. 

It can be argued that constraining the dark matter interaction based on a single spuri- 
ously small black hole is unreasonable considering other effects that could conspire to produce 
a small black hole in the scenario described in §2. For example, the black hole in M33 could 
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have been ejected in a merger event, or more importantly, the value of a in M33's dark halo 
could have been small. Figure 1 indicates that the black hole mass is extremely sensitive 
to a for a given physical cross section for a = — 1. Since the real cosmic variance of a is 
certainly at least 0.1 — 0.2 (Subramanian et al. 1999), it is plausible that a smaller than av- 
erage value of a in the post-collapse dark halo of M33 is responsible for the small black hole 
in this system. However, M33 is the best observed local example of a group of several bulge- 
less galaxies that lack supermassive black holes comparable to those found in bulge systems 
(Richstone et al. 1998; Kormendy & Gebhardt 2001): NGC 4395 (Sm) M BH < 8 x 10 4 M Q 
(Fillipenko & Ho 2001), IC 342 (Scd) M BH < 5 x 10 5 M (Boker et al. 1999), NGC 205 
(dE5) M BH < 9 x 10 4 M Q (Jones et al. 1996). It is doubtful that cosmic variance in a can 
account for the small black holes in all of these systems. 



5.2. Black Hole Constraint from Bulge Systems 

An even tighter constraint on the dark matter interaction can be obtained from bulge 
systems if compression of the dark halo is allowed for. In Appendix B we show that conden- 
sation of baryons results in a steeper inner density profile with a' = /(«,£), where a and a' 
are the inner dark matter profile exponents before and after compression respectively, and 
£ is the exponent of the inner profile of the total mass density (baryons plus dark matter) 
after compression. Specifically, we show that, in both the collisionless (A/r ^> 1) and highly 
collisional (A/r <C 1) limits, the final dark matter cusp slope is 

a=f{a,0 = - A • (5-2a) 

4 — a 

For a final flat rotation curve, typical of normal bulge dominated galaxies, £ = 2, giving the 
simple result 

/ 6 — a . _ . 

a = . 5-2b 

4 - a 

The most likely CDM value of a = 1.3, gives a' = 1.74. The tightest constraint is obtained 
by comparing to the smallest bulge system that lies on the observed M BH — v relation, since 
roughly speaking the M BH — v scaling for SIDM black holes is as steep as or shallower than 
that observed. To this end, we compare to the black hole at the galactic center, M BB ~ 
3.0 x 10 6 M© (Genzel et al. 2000), which is the smallest reliable mass estimate that lies on 
the observed M BB —v (Ferrarese & Merritt 2000; Gebhardt et al. 2000; Merritt & Ferrarese 
2001). Using this value in eqn. (5-1), with r s and p s determined from the scaling relations 
for Mhato = 3.0 x 10 12 M gives, 

a lV a W0 < 0.02 (0.92) a (5-3) 
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for a' = 1.74. Our "best fit" model (a = 1.7, a = 0, and <Tif" 00 = 0.02), mentioned in 
§4, will obviously satisfy the constraint in eqn. (5-3), since it was chosen to reproduce the 
magnitude and slope of the observed M BH — v relation. 

6. Mergers 

In the black hole formation scenario described in this paper, we have considered the 
growth of a supermassive black hole from dark matter at the center of a single dark halo 
from a stellar mass size seed. However, in standard hierarchical structure formation models, 
massive halos experience multiple mergers during their lifetimes. This suggests that the 
mass of the supermassive black hole in a galaxy at z = could be the sum of smaller black 
holes grown in its progenitor halos. From eqn. (2-8), an efficiency for growing black holes 
in a single dark halo, e (Mhaio) = Mg ^^ ai °' > , can be calculated. For the region of parameter 
space which satisfies all constraints (see below), < a < 1, we find that e oc M^ alo , where 
—0.1 < p < 1.4 for a 1.3 (—0.1 < p < 0.5 for a 1.7), and larger values of p correspond 
to smaller values of a. Thus, the efficiency is roughly a constant for a = 1, and varies as 
a small positive power of halo mass for a < 1. While a detailed calculation of halo merger 
dynamics is beyond the scope of this work, it should be noted that because e (Mhaio) is a 
weak function of Mhaio for the models considered, we make little error in simply using the 
present dark halo to calculate the black hole mass. 

The fact that our scenario populates dark halos at high redshift with supermassive black 
holes has other interesting astrophysical consequences. Consider our "best fit" model with 
a = 1.7, a = 0, and <7iVi 00 = 0.02. For this set of parameters the efficiency for growing black 
holes in a given dark halo is approximately a constant, e cx M^ a , so that the black hole in 
each dark halo will be proportional to the mass of the halo. The collisional cross section for 
this model is so small that collisional relaxation will not ameliorate the cuspy halo problem, 
as the relaxation timescale will be much longer than the Hubble time. However, N-body 
simulations of mergers of galaxies containing supermassive black holes have demonstrated 
that mergers between galaxies with steep power-law density cusps produce remnants with 
shallower power-law cusps, because the formation of a black-hole binary transfers energy to 
the halo, lowering the central density (Merritt & Cruz 2001; Milosavljevic & Merritt 2001). 
Consider then, the early generations of progenitor dark halos at high redshift that grew 
black holes by accretion of SIDM. This generation might correspond to the epochs at which 
the first massive stars collapsed to form seed black holes (Madau & Rees 2001; Abel et al. 
2000, 2001). In accordance with our calculation in Appendix B, we expect an inner profile 
cusp of a ~ 1.7, in these systems, and black holes will be grown in each dark halo with 
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an efficiency e. Successive mergers between progenitor halos harboring supermassive black 
holes will result in a milder density cusp in each remnant, so that parent halos at z = will 
have significantly shallower density profiles, thus providing a possible solution to the cuspy 
halo problem. 

As inner profile cusps become shallower at late times, there will come a point at which 
black holes only grow by merging: accretion of SIDM will be negligible because of the steep 
dependence of Mbh on a. Furthermore, because the efficiency, e, is roughly constant, the 
black hole masses will scale in proportion to the masses of their host halos. Then, as pointed 

12 

out by Haehnelt and Kauffmann (2000), the scaling Mhaio oc v s (n_1) for power law power 
spectra (see §4), gives Mbh oc at galaxy scales where n R3 —2. Given this scaling initially, 
mergers between galaxies will move black holes along the observed Mbr—v relation. Finally, 
it should be noted that the dark matter cusp after compression a', will have much less scatter 
than the initial profile, since a range of 1 < a < 2 is compressed to 1.67 < a' < 2. The 
corresponding black hole mass will also be less sensitive to a' as the curves in Figure 1 are 
flatter for larger a. It follows that the efficiency e and M B h — v scaling are determined by a 
set of fixed parameters: the fundamental physics of the dark matter interaction determines 
Ci^ioo an d ct and the universal (initial) profile of dark halos with a = 1.3 ±0.2 coupled with 
the adiabatic collapse of baryons (see Appendix B) determines a' ~ 1.7 — 1.8. 

Accordingly, the black hole growth scenario with the parameters of our "best fit" model 
has the following three consequences: 1) The observed magnitude and scaling of the observed 
Mbh ~ v relation (at z = 0) is reproduced. 2) The lack of comparable supermassive black 
holes in bulgeless galaxies like M33 is explained 3) It provides a possible solution to the 
cuspy halo problem via mergers. 

This scenario becomes all the more plausible when one realizes that massive black holes 
need only grow in a few percent of progenitor dark halos. Menou, Haiman, & Narayanan 
(2001), have shown that the presence of central massive black holes at the centers of nearly all 
nearby galaxies can arise from their merger history, even if only a small fraction (~3x 10~ 2 ) 
of the progenitor halos harbored black holes at high redshift. Hence, even if there were 
significant scatter in the final inner profile slope a', which would yield a corresponding 
scatter in black hole mass (because of the steep dependence of M B h on a) , our model would 
still have the aforementioned consequences as a result of the merger history of dark halos. 
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7. Constraints from Other Observations 

The evolution implied by the collisional nature of dark matter must not conflict with 
observations of dark halos over nearly three orders of magnitude in halo circular velocity. In 
this section we discuss four constraints on the dark matter interaction obtained at different 
halo mass scales. First, a minimum cross section per unit mass is required to solve the cuspy 
halo problem on dwarf galaxy and low surface brightness (LSB) galaxy scales, below which 
SIDM interpolates smoothly with CDM and is astrophysically uninteresting. Second, all dark 
halos observed today must not have undergone core collapse by the present, giving a upper 
limit or core collapse constraint. Third, as pointed out by Gnedin and Ostriker (2001), 
an upper limit on the dark matter interaction can be obtained by requiring that galactic 
subhalos survive until the present in hotter cluster environments. Finally, an observational 
upper limit on the core radius of a cluster of galaxies (Arabadjis et al. 2001) also gives a 
strong constraint on the dark matter interaction. 



7.1. Uninteresting Limit and Core Collapse Constraint 

In the regime where SIDM halos are optically thin, the dynamics resembles two body 
relaxation in globular clusters. Heat will flow inward from the outer halo due to the tem- 
perature inversion implied by the post-collapse profile (see Figure 2). A flattened core will 
develop and grow outward, with the central density falling and the velocity rising. Once the 
temperature inversion is gone, expansion halts. The direction of heat flow reverses and core 
collapse begins. The state of evolution of a dark halo will be determined by the ratio of the 
relaxation time at its characteristic scale to its age. For example, if core collapse begins after 
a number C\ relaxation times, then all halos observed today must have 

C, t rel (r s ) = ° l Uyn {rs) >t H - t f , (7-la) 

where r s is given by eqn. (2- 7b), t# is the age of the universe, and tf is the halo formation 
time. Furthermore, if a number C2 relaxation times must pass before the central density 
cusp is significantly flattened, then SIDM will not resolve the cuspy halo problem unless 

C 2 t rel (r s ) <t H -t f . (7- lb) 

Or, 

^<^<^, (7-lc) 

a constraint which must be satisfied by all non-cuspy halos. The constants C\ and C2 can 
be calibrated against simulations. 
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The formation times in eqn. (7-1) are calculated by taking the median of the formation 
time distribution of Lacey & Cole (1993), which is based on the extended Press-Schecter 
formalism (Press & Schecter 1974; Bond et al. 1991). The formation time of a halo of 
present mass M is defined as the time when a parent halo appeared which had half or more 
its mass. This particular definition of formation somewhat remedies the neglect of accretion 
and merging in eqns. (7-la) and (7- lb), since the countdown to core collapse starts only 
after the last major merger in a halos formation history. Table 1 lists formation times and 
redshifts of dark halos from the dwarf to cluster scale. These have been calculated using the 
analytic fitting formulae in Appendix B & C of Kitayama & Suto (1996) for the CDM power 
spectrum (Bardeen et al. 1986) and formation time distribution (Lacey & Cole 1993). 

Note that the relaxation time and age of a halo will depend on the scale under consider- 
ation, so that systems at different scales will be at different stages of evolution. Specifically, 
the velocity dependence of the cross section in eqn. (2-6) implies lower collision rates in 
larger systems, so that, for example, dwarf galaxies will be more relaxed than clusters. Fig- 
ure 4 plots the inequalities in eqns. (7-la) and (7-lb) in the £7ii>i 00 - v C i rc plane, where v C i rc 
is the circular velocity of the dark halo, for several values of a and values of C\ and C2 
taken from simulations. The leftmost and rightmost curves correspond to eqn. (7-lb) and 
(7-la) respectively. Each set of parallel curves corresponds to a different value of a and the 
region between each set of curves satisfies both inequalities. To the left of the leftmost curve 
inequality (7-lb) is violated: <Jivf 00 is too small and SIDM halos are indistinguishable from 
CDM halos. To the right of the rightmost inequality (7-la) is violated: <7ii>" 00 is too large 
and dark halos have core collapsed. At fixed v circ the region between the curves indicates the 
range of (Jivf 00 for which halos at that scale will be significantly evolved due to collisions. Or 
for fixed (JiVi 0Q , it indicates the range of scales that will have undergone the desired amount 
of evolution. 

From the curves in Figure 4 for a = 0, a clear but somewhat surprising result emerges: 
the quantity t ^ r l l t is nearly a constant over three orders of magnitude in circular velocity. 
A constant cross section thus implies dwarf galaxies will be just as relaxed as clusters. 
However, observations indicate that dwarf galaxies have larger cores in proportion to their 
characteristic size than do clusters. There is thus a scaling problem for a constant cross 
section as pointed out by Dave et al. (2000) and Yoshida et al. (2000b). A velocity dependent 
cross section obviously remedies this problem since collisions will be more frequent in lower 
velocity environments and smaller systems will thus be more evolved. However it is apparent 
from Figure 4 that if the velocity dependence is too steep (a is too large) only a narrow 
window of scales will have evolved significantly and all larger scales will be identical to CDM 
halos. We can obtain a rough estimate of the value of a that will reproduce observations of 
core sizes on both dwarf galaxy and cluster scales as follows. Kochanek & White's (2000) 
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simulations indicate that the core radius grows linearly during the expansion phase of halo 



Consider DDO 154, a dwarf galaxy which has a core radius r c ^ w ~ 2.5 kpc and maximum 
rotational velocity v max ~ 47 km s _1 (Carnigan & Purton 1998), which can be identified 
with the circular velocity of the dark halo. The lensing cluster EMSS 1358+6245 has a core 
radius r CjCi < 40 kpc and an inferred halo mass of M ha i Q 4 x 10 14 M Q (Arabadjis, Bautz, & 
Garmire). Using the scaling relations in Appendix A to calculate p s , r s , and v s , eqn. (7-2) 
gives a w 0.6. 

The curves in Figure 4 can be translated into constraints on the dark matter interaction 
by fixing v circ and thus specifying a physical scale. For example, since smaller mass halos 
will always be more evolved relative to larger halos, an upper limit can be obtained in 
the (JiV^qq — a plane by requiring that the smallest dark halos observed today have yet 
to undergo core collapse. We designate v C i rc = 20 km s _1 halos as the smallest observed 
today and generate a core collapse constraint which is plotted in Figure 5. Similarly, a 
lower limit on the dark matter interaction is obtained by requiring that v c i rc = 50 km s _1 
halos, corresponding to dwarf and LSB galaxies, have undergone significant evolution by 
the present. Below this limit, shown in Figure 5, SIDM effectively becomes CDM and is 
astrophysically uninteresting. 

These constraints depend upon the constants C\ and C<i in eqns. (7- la) and (7- lb) 
respectively, which, strictly speaking could be functions of a since the dynamics could differ 
slightly for different cross section velocity dependences. However, accurate simulations have 
only been carried out for a constant cross section so we take them as constants and calibrate 
to simulations for a = 0. Quinlan's (1996) Fokker-Planck simulations indicate that density 
profiles containing temperature inversions enter a core collapse phase after approximately 

5 half- mass relaxation times. Burkert's (2000) N-body simulations roughly agree with this 
result, but Kochanek & White (2000) saw much faster evolution towards collapse. Balberg 
et al. (2001) showed that the core collapse time of an SIDM halo is 290 central relaxation 
times, however this is based on the assumption that the halos have flat inner profiles at 
all times, which, as per our discussion at the beginning of §2, is highly unlikely to hold. 
We took Ci = 4.76 in agreement with the core collapse timescale determined by Kochanek 

6 White (2000), converted to our units. Simulations indicate that a minimum value of 
cdm = 0.45 cm 2 g" 1 is required to flatten dwarf galaxy and LSB galaxy halos (Dave et 
al. 2000; Wandelt et al. 2000). Plugging this cross section into eqn. (7- lb) at the scale 




. Then the ratio of cluster to dwarf core radii, r^ c [/r^ dw , 



Log 



tH-tf,d (psrsVs) cl r Ctdw 
tH~tf,dw (Psr s v s ) dw r c>ci 



(7-2) 



a = 



Log {v StC i/v s4w ) 



-19- 



Vdrc = 50 km s 1 with a = gives C 2 = 0.48. 

7.2. Evaporation Constraint 

The existence of subhalos in larger halos further constrains the dark matter interaction. 
For galactic subhalos in clusters, heat transfer to a cool subhalo from the hot cluster en- 
vironment could evaporate the galactic halo if heat transfer were too efficient. This would 
violate the fundamental plane relations in conflict with observations. Gnedin and Ostriker 
(2001) excluded a range of dark matter cross section per unit mass by requiring that galactic 
subhalos survive until the present epoch. They considered both optically thick and thin halos 
and obtained lower and upper limits respectively on the dark matter interaction strength. 
Both upper and lower limits can be obtained because of the non-monotonic dependence of 
the heat transfer on the interaction cross section. We concern ourselves only with the upper 
limit obtained in the optically thin regime, since, as mentioned previously cross sections large 
enough to give optically thick halos can be excluded on other grounds (Dave et al. 2000; 
Wandelt et al. 2000). To obtain their scattering regime constraint, Gnedin and Ostriker 
required the cluster relaxation time evaluated at the position of the galactic halo be longer 
than the age of the cluster, which they approximate as the Hubble time. A typical galactic 
halo will be located at the scale radius of the cluster r s , so that the evaporation constraint 
is just given by eqn. (7-la). As a conservative estimate, we take the constant C± — 2 for 
evaporation. This constraint is plotted in Figure 5 for a cluster like EMSS 1358+6245 with 
a halo mass M halo «4x 10 14 M© (Arabdjis et al. 2001). 

7.3. Cluster Core Constraint: Cluster EMSS 1358+6245 

Recently Arabadjis, Bautz and Garmire (2001) used high resolution Chandra obser- 
vations of the lensing cluster EMSS 1358+6245 to constrain the dark matter interaction 
strength. Specifically, an upper limit of 42 kpc on the size of any constant density core in 
the M ha i ~ 4 x 10 14 M cluster was obtained from their X-ray determination of the cluster 
mass profile. This upper limit was then compared to Yoshida et al.'s (2000b) numerical 
simulation of cluster size SIDM halo, where a comparable mass cluster had a core ~ 40 kpc 
for (T = 0.1 (a = 0). This allowed Arabadjis et al. (2001) to place an upper limit of er < 0.1 
(a = 0) on the scattering cross section. In the context of the discussion in this section, 
Arabadjis et al.'s constraint is equivalent to evaluating eqn. (7-la) for M^aio ~4x 10 14 M , 
o"o — 0.1, and a = 0, thus determining the number of relaxation times (Ci) it takes for 
the cluster to develop a ~ 40 kpc core. This equation is used to plot the lensing cluster 
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constraint in Figure 5. 

8. Conclusion 

The shaded area in Figure 5 indicates the region of the (JiVi 00 — a plane consistent with 
the four constraints discussed in the previous section. A window of possible dark matter 
interactions with 0.5 < a < 3 and 0.5 < Ci^ioo ~ 5 appears to satisfy the uninteresting, core 
collapse, evaporation, and cluster EMSS 1358+6245 constraints, which is broadly consistent 
with the results of numerical simulations (for the regions of parameter space that have been 
simulated). 

In Figure 6 the permitted region of parameter space is plotted with the black hole 
constraints from eqns. (5-1) and (5-3) for both bulgeless (M33) and bulge (Milky Way) 
systems. The bulgeless constraint is plotted for three plausible values of the inner profile 
slope (without baryonic compression) a = 1.1, 1.3, and 1.5 from right to left. The region 
to the left of each respective curve is permitted by the constraint. The bulge constraint is 
plotted for a final profile cusp of a' = 1.74 as determined in Appendix B for an initial profile 
cusp with a most probable value of a = 1.3 (even if a = 1.0, the NFW value, the final cusp 
is a' = 1.67). The region to the left of the nearly vertical curve is permitted. 

From the bulgeless constraint in Figure 6, it is apparent that if the post-collapse density 
profile of the dark halo of M33 had a > 1.3, no region of the parameter space of dark matter 
interactions exists that is consistent with both the black hole in M33 and constraints from 
other observations. This would exclude SIDM as a possible solution to the problems with 
CDM on subgalactic scales. If a = 1.3, a tiny region is still allowed with a « 0.5 and a x v " 00 « 
0.5, corresponding to lower vertex of the triangle in Figure 6. A much larger region of 
parameter space would be available for a < 1.1. The gray shading in the figure indicates the 
region of parameter space allowed if a — 1.1 for M33; however, N-body simulations of dark 
matter clustering indicate that cusps this mild are unlikely. Apparently, the cross section 
velocity dependence a = 1 discussed extensively in the literature can be excluded. Only 
detailed N-body simulations can determine whether a dark matter cross section consistent 
with these regions of parameter space can reproduce the observed structure of dark halos 
on all scales. With regards to black holes, the parameters quoted above with a = 1.3, 
a « 0.5, o"it>i 00 ~ 0.5, would account for < 1% of the supermassive black hole mass at the 
centers of galaxies. This is roughly the ratio of the observed upper limit on the black hole 
in M33, M BH < 1500 M , to the M BH ~ 5 x 10 5 M that M33 would have if it lay on the 
observed Mbh — v relation. This correspondence arises of course because we have forced 
these parameters to be consistent with the black hole constraint from M33. Although this 
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fraction is small, the resulting black hole masses (typically > 10 4 for normal galaxies) are 
far more massive than stellar mass black holes, and are appropriate seeds for rapid baryonic 
growth. 

The constraint from bulge systems, which takes into account the probable compression of 
the dark halo by baryons, is much tighter than the bulgeless constraint, effectively requiring 
<J\v1 m < 0.02 for a < 1. Interaction strengths consistent with this constraint are too small 
to produce significant collisional evolution of dark halos, since t re i ^ — tf, so that in the 
absence of other dynamical processes, SIDM behaves effectively as CDM, failing to remedy 
CDM's small scale problems. Nevertheless, the interaction parameters in the range permitted 
by the bulge constraint, would still seed dark halos with supermassive black holes at high 
redshift. As mentioned in §6, SIDM could then provide an indirect solution to the dark halo 
cusp problem, since successive mergers of dark halos harboring supermassive black holes at 
high redshift would result in remnants with flatter inner density profiles. 

Balberg & Shapiro (2001) have considered an alternative supermassive black hole for- 
mation scenario, whereby supermassive black holes are formed after SIDM halos undergo 
gravothermal collapse. However, for the range of dark matter cross sections permitted by 
the constraints in Figure 6, our calculations of the core collaspe timescale in §7 and depicted 
in Figure 4, indicate that core collapse is highly unlikely for dark halos with circular velocities 
between 10-1000 km/s. 

In sum, the self interacting dark matter scenario remains of great interest. If the dark 
matter interaction is a steep function of velocity (a > 1), a broad range of parameter space 
with 1.0 < a < 3 and 0.5 < o"ifi 00 ^ 5, is available. Though, as mentioned previously, 
a steep velocity dependence necessarily implies only a small range of dark halos masses 
will show significant collisional evolution (see Figure 4). For shallower velocity dependences 
(a < 1), the range of astrophysically permitted interaction strengths imply heat conduction 
driven flattening of dark halo profiles is unimportant. However, the early and efficient growth 
of massive black holes at the centers of dark halos is a natural consequence of the theory, 
which in turn, through succesive mergers, leads to the observed density profiles in galaxies, 
produces black holes with the appropriate mass and mass scaling, and explains the lack of 
massive black holes at the centers of bulgeless galaxies like M33. 
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Martin Rees, Uros Seljak, David Spergel, Paul Steinhardt, and Scott Tremaine for helpful 
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A. ENS SCALING RELATIONS 

The post collapse density profile of an SIDM dark halo will be identical to that of a 
CDM halos. High resolution numerical simulations (Navarro et al. 1996,1997) have shown 
that these profiles may be fitted to a universal shape, which we take to be a Zhao profile 
(Zhao 1996) 

P(r) " ) \ 6 — ol ' (Al) 

\r s J (l + r/r s ) 

Here r s is a characteristic length scale and p s = 5 C p cr u is a characteristic density, which is 
equal to a density enhancement 8 C times the critical density for closure p cr u = 3H 2 /8nG. 
The two free parameters 5 C and r s can be determined from the halo concentration ca (Ma) 
and the virial mass Ma 

A r 3 

$c = - —, (A2) 

3 Zra(l + c A ) -c A /(l + c A ) V ; 

r A 9.52 x 1(T 2 f M- x 1/3 



ca ca 



ij£) ^ 2/3A " /3k P c (A3) 



where 



qo.30 if A = 

a (n , a) = 178^; i A n + l = i- (A4 > 



Here is the virial radius, inside which the average overdensity is A times the critical 
density for closure, and Ma is the mass within For any particular cosmology, the 
concentration ca (Ma) is a function of the virial mass which results from the fact that dark 
halo densities reflect the density of the universe at their formation epoch, and smaller mass 
halos collapse earlier in hierarchical structure formation. Hence, for CDM power spectra, 
ca (Ma) is a decreasing function of Ma- ENS (2001) have carried out an extensive suite of N- 
body simulations to characterize the dependence of ca (Ma) on the cosmological parameters 
by fitting their simulated dark halos to an NFW profile, which is eqn. (Al) with a = 1 and 
e = 3. We have used their publicly available routine to calculate ca (Ma) for the ACDM 
cosmology quoted in §1. Strictly speaking, these concentrations apply only to the NFW 
profile; however, we use them to characterize the more general profile in eqn. (Al), since 
ca (Ma) should not change significantly if the inner profile slope a is instead taken as a free 
parameter. 
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Given the characteristic radius and density, r s and p s , the characteristic velocity dis- 
persion v s can be calculated from v 2 s = fiGp s r 2 where \i = / 3 _au a -i) ■ The density profile 
in eqn. (Al), will give a circular velocity profile, v C i rc (r), which increases outward from the 
center passing through a maximum at r ~ r s and then decreases. The circular velocity of 
a dark halo is conventionally taken to be the value at this maximum, since r s is the only 
length scale in the system. For the NFW profile, this maximum occurs at r th 2r s , but for 
our more general profile the location of the maximum will depend on a and e. We define the 
circular velocity of a dark halo v circ , to be the maximum of the NFW profile (eqn. (Al) with 
a = 1.0 and e = 3.0), again assuming that v C i rc will not change significantly if a is allowed 
to vary. 

v circ = 3.13 x 10- 3 ( ( — / A ^ V2 h-^A 1 ^ km s" 1 (A5) 

\M e J \ln(l + c A )-c A /(l + c A )J 1 ; 



B. DENSITY CUSP ENHANCEMENT BY DISSIPATIVE COLLAPSE OF 

BARYONS 

The dissipative infall of baryonic matter will strongly perturb the underlying dark matter 
distribution, pulling it inward and resulting in a steeper density cusp in the innermost 
regions of the dark halo. If the dark matter is collisionless and the baryons condense slowly 
(compared to the dynamical time in the inner halo), the angular momentum, or equivalently 
rM (r), is an adiabatic invariant, allowing one to calculate the final density profile of the dark 
halo given the initial dark matter profile and the final profile of the baryons (Blumenthal et 
al. 1986; Flores et al. 1993). 

For our purposes, we are interested in the density profile of the dark halo for < r < r*, 
where r* is defined as the radius of optical depth unity: 

n=r s Tf^, (Bl) 

with the optical depth given by eqn. (2-7a) and r s by eqn. (2-7b). For the density profiles 
given by eqn. (2-1), the relaxation time goes to zero at the center. This sets a limit 
to the minimum radius, rn, for which the system can be considered collisionless, where the 
relaxation time equals the age of the halo. For the range of scales and parameters we consider, 
r* <C th] hence, dark matter is highly collisional in the region where the density profile is 
desired, and the adiabatic approximation breaks down: collisions will scatter dark matter 
particles onto orbits of different angular momentum, breaking the invariance of rM (r). 
For regions of the halo with r (r) ~ 1 (r ~ r*) the adiabatic approximation should hold 



reasonably well; however, for regions with r (r) 3> 1 (r <C r*), the dark matter behaves as a 
fluid, so that the equations of hydrostatics must be employed. 

In what follows we determine the steepening of the inner density cusp in the collisionless 
and highly collisional (fluid) limits, and find the slope of the final dark matter density cusp to 
be identical. We do not know of a simple approximation to determine the cusp in the region 
where the dark matter is moderately collisional, but expect the density profile to interpolate 
smoothly between the inner (fluid) and outer (collisionless) regions of the halo, as the cusp 
slope will be the same. In the following derivations we drop all numerical constants and 
consider only the scaling of various quantities with r. Subscripts % and / are used where 
appropriate to distinguish the initial and final states. In the initial state, the dark matter 
and baryons are assumed to be well mixed, so that the density profile of the baryons will 
be parallel to that of the dark matter. In the final state, the baryons are assumed to 
dominate the total mass (dark matter plus baryons) Mf (r/) in the inner region where the 
dark matter profile is desired, and this final total mass distribution is taken to be Mf ~ 
corresponding to a total density profile pf ~ rj^. We assume the dark matter density cusp 
is pi ~ r~ a before condensation and pf ~ rj a after, and determine a' = f (a, £). 

For the collisionless case, we make the simplifying approximation that the dark matter 
particles are on circular orbits, which is justified since there is more phase space available 
for nearly circular orbits than for radial ones (Blumenthal et al. 1986; Flores et al. 1993; 
Navarro, Frenk & White 1996). Conservation of dark matter implies that 

p^dn = p f r 2 f dr f rf~ a ~ rj- a ' . (B2) 

Conservation of angular momentum gives 

TiMi (ri) = TfMf (rf) ~ =>- rf~ a ~ r 4 ^^ 

Combining these two results we find 

, a + 3£ - a£ 

a = a 

4 - a 



(B3) 



(B4) 



In the fluid limit, we assume the dark matter is in hydrostatic equilibrium in the initial 
and final states. We neglect heat transfer, so that the entropy of the dark matter at fixed 
mass shell is invariant 

Si {n) = S f (r f ) , (B5) 

where S is given by 

P(r) 
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for a monatomic ideal gas, and r\ and r/ are related by mass conservation (eqn. (B2). These 
assumptions are valid provided the condensation of baryons takes place on a timescale slow 
relative to the dynamical time, but fast compared to the heat conduction timescale (eqn. 
(3-4)). If ts is the condensation time, these approximations hold where td yn (r) < ts < 
t ( r ) tdyn (r). Hydrostatic equilibrium of the dark matter component implies 



1 dP M (r) 



P ~ pMr- 1 ; (B7) 



p dr r 2 
or 

Pinrf-**; P f ~r)-Z- a '. (B8) 
Combining with eqns. (B5) and (B6) gives 

2-a/3 2-£+2a'/3 / Dn \ 

r { ~?/ , (B9) 

and mass conservation (eqn. (B2) then implies 

ol = + (BIO) 
4 - a 

which is identical to the collisionless case. 

For the most probable value of a — 1.3, and taking £ = 2 for a final flat rotation curve, 
typical of normal bulge dominated galaxies, a' = 1.74. Note that in deriving our black 
hole constraint from bulge systems in §5 we use the initial state values of r s , p s , and v s , 
since compression will only change these quantities by factors of order unity, resulting in an 
insignificant change in Mbh compared to the change of several orders of magnitude caused 
by the steeper cusp a' (see Figure 1). 
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TABLE 1 

Formation Redshifts and Times for Dark Halos 



M halo 


z f 


tf/t H 


(M e ) 






10 8 


1.83 


0.26 


10 9 


1.61 


0.29 


10 10 


1.38 


0.33 


10 11 


1.15 


0.39 


10 12 


0.93 


0.45 


10 13 


0.72 


0.52 


10 14 


0.53 


0.61 


10 15 


0.36 


0.71 



NOTE.— Formation times and redshifts for dark halos from the dwarf galaxy to cluster scale 
in a ACDM cosmology. These are calculated from the Lacey & Cole (1993) formation time 
distribution. M ha i Q is the virial mass of the dark halo, Zf the formation redshift, and tf/t H 
the ratio of the formation time to the age of the universe. 
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Fig. 1. — Black hole mass grown from optically thick accretion as a function of the inner 
density profile exponent a for different values of a (the velocity dependence of the cross 
section) for a Milky Way sized dark halo. Solid, dotted, and short-dashed curves are for 
a = 0,0.5, 1 respectively. The dark matter interaction strength has been set to cif" 00 = 1. 
The Mbh is extremely sensitive to a for a — — 1. 
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v^(r)=3k B T/2m p 




Fig. 2. — Schematic illustration of the density and velocity dispersion (temperature) profiles 
for SIDM dark halos at early times. A black hole accretes dark matter from the center of 
the halo. Power law slopes are indicated near each segment of the profiles. Arrows indicate 
the direction of heat trasfer made possible by the collisional nature of the dark matter. 
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Fig. 3. — Black hole mass grown from optically thick accretion as a function of dark halo 
circular velocity for different values of a (the velocity dependence of the cross section) for 
a Milky Way sized dark halo. The shallow and steep dotted lines are the observed scaling 
from Gebhardt et al. (2000) and Merritt & Ferrarese (2000) respectively. The solid lines are 
the M B h — v relation for SIDM with a = 0, 0.5 and 1, from steepest to shallowest (a = 1 is 
the shallowest curve). The dark matter interaction strength has been set to o"ifi 00 — 1 an d 
the inner profile exponent is a = 1.5. The dashed curve is our "best fit" model (a = 1.74, 
a = 0, <7iu?oo = 0.02) 
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Fig. 4. — Plot of the inequalities in eqns. (7-la) and (7-lb) for C2 = 0.48 and C\ = 4.76. The 
leftmost and rightmost curves correspond to eqn. (7-la) and (7-lb) respectively. Each set of 
parallel curves corresponds a different value of a and the region between each pair satisfies 
both inequalities. Solid, dotted, and dashed lines refer to a = 0, 1, and 3 respectively. 
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Fig. 5. — Constraints on the dark matter interaction plotted in the civ® 00 — a plane. The 
solid line is the uninteresting limit, corresponding to the minimum cross section required to 
flatten dwarf and LSB galaxies. The dashed curve is the core collapse constraint, obtained 
by requiring that 20 km s -1 dark halos have yet to undergo core collapse. The dotted line is 
the upper limit adapted from the evaporation constraint of Gnedin & Ostriker (2001). The 
constraint obtained from Chandra observations of cluster EMSS 1358+6245 (Arabdjis et al. 
2001) is the short dashed line. Hash marks indicate the region excluded by each constraint. 
The shaded region is consistent with all four constraints. 
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Fig. 6. — Black hole constraints on the dark matter interaction plotted with the available 
region of parameter space for SIDM. Dotted, solid, and short dashed lines are the black hole 
constraints obtained by comparison to the bulgeless galaxy M33 for a = 1.1, 1.3, and 1.5 
respectively. Hash marks indicate the region excluded by each constraint. The large triangle 
is the region of parameter space allowed by the constraints discussed in §6 and shown in 
Figure 5. Shading indicates the region which is consistent with the constraints from §6 and 
the M33 constraint with a = 1.1. The long dashed curve is the black hole constraint obtained 
by comparing to the Milky Way after taking into account the adiabatic compression of the 
dark halo by baryons, which changes the inner profile slope to a' = 1.74. The region to the 
left of this curve is permitted. 



